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ABSTRACT 

We present new observations of two anomalous X-ray pulsars (AXP) 4U 0142+61 
and IE 1048.1—5937 made in 1998 with the ASCA. The energy spectra of these two 
AXPs are found to consist of two components, a power-law and a blackbody emission 
from the neutron star surface. These observations, when compared to earlier ASCA 
observations in 1994 show remarkable stability in the intensity, spectral shape and pulse 
profile. However, we find that the spin-down rate in IE 1048.1—5937 is not constant. 
In this source, we have clearly identified three epochs with spin-down rates different 
from each other and the average value. This has very strong implications for the mag- 
netar hypothesis of AXPs. We also note that the spin-down rate and its variations in 
IE 1048.1—5937 are much larger than what can normally be produced by an accretion 
disk with very low mass accretion rate corresponding to its low X-ray luminosity. 

Subject headings: stars: neutron — Pulsars: individual (4U 0142+61, IE 1048.1—5937) 
— X-rays: stars 
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Introduction 



Some X-ray pulsars are known lo have remarkable 
similarity in their properties which are dif ferent from 
other binary or iso lated X-ray pulsars ( Mereghetti 
'k and Stella 1995). The properties common to 



mo et of these objects arc a) pulse period in a small 



ran ge of 5 12 s, b) monotonous spin down with P/P 
in the range of 5 x 10^^ — 1.3 x 10^'^ s c) iden- 
tical X-ray spectrum consisting of steep power-law 
(r = 3 — 4) and black body component (kT 0.5 
keV), d) stable X-ray luminosity (10^'' — 10'^^ ergs 
s~^) for years, e) faint or unidentified optical coun- 
terpart, and f) no evidence of orbital motion. The 
sources also have a galactic distribution, most of 
these are within |6| < 0.5° and all are probably 
young (~ 10^ yr) because of their association with 
SNR or molecular clouds. The objects in which all 
the properties mentioned above have been observed 
are 4U 0142-^61, IE 2259+586, IE 1048.1-5937, 
IRXS J170849.0-400910 and IE 1841-0 45 (Kes 73). 
Two more objects, AX J1845.0-0300 ([Torii et al 



19S|8 ) and RX J0720.4-3125 (Ilabcil el al. 1990 ^ 



als o probably belong to the same class but to estab- 
lish their AXP candidacy, more X-ray observations 
are required to measure their pulse period variations, 
search for possible pulse arrival time delay and inves- 
tigate the fiux stability. Classification of an object as 
AXP only from some properties similar to the above is 
not very firm. 4U 1626—67, probably a binary system, 
showed both spin- up and spin-down (Chakrabarty et 
il. 1997) and also has an opt ically bright accretion 



disk ( Middleditch et al. 1981). Ano ther one object 



RX J1838.4-0301 ( |Schwentker 1994| ), does not have 
stable intensity and pulsations are also not always 
detectable ( ^ong et al. 199£ ). Therefore, these two 
objects are not AXPs. 

Considering the strong similarity between these 
ha ndful of sources, it has been proposed that they 



ha ve same physical nature and different aeenarioa 
have been proposed to explain the observed proper- 
ties. The prominent models are a) accretion from 



lowi 



biuai)' Lumpauiun (MeiegheLLi &l and GLelld 



19Ste), b) single ueuLiuu sLai acLiuLiug hum uiulecu- 



lar cloud, or a product of common envelope evolution 
(Thorne-Zytkov object) of close high mass X-ray bi- 
naries in which a solitary neutron star accretes mat- 
ter from a fossil disk (van Paradijs, Taam, & van den 
Heuvel [1995| ; Ghosh, Angehni, & White |1997D , and 
c) extremely high magnetic field neutron star radi- 



ating X-rays due to magnetic field decay (Thompson 
& Duncan 1996). Unlike the radio pulsars and ro- 



tationally powered X-ray pulsars, in the AXPs, the 
spin-down rate is not large enough to power the ob- 
served X-ray emission. 



Among the 90 or so known X-ray pulsars (Nagase 



1999), direct evidence of binary nature is known for 



more than 35 sources (van Paradijs 1995). Including 
the indirect evidences this number can be upto about 
65 and 7 pulsars are isolated stars in SNR and are 
powered by rotational energy losses. In the rest of 
the pulsars, in which no binary signature is known, 
it is often due to lack of sufficient observation. How- 
ever, the 7 objects which are either AXPs or candidate 
AXPs (or 10 if we include the 3 Soft Gamma-ray Re- 
peaters in which pulsations have been detected), no 
binary signature has been found in spite of extensive 
searches. The strong upper limit on pulse arrival time 
delay that has been obtained in some of these sources 
strongly suggests non-binary nature for the AXPs. In 
addition, the spin change behaviour of the AXPs is 



also remarkably different from accreting pulsars (Bild 



sten et al. 1997). In almost all the accreting X-ray 



pulsars, both spin-up and spin-down episodes have 
been seen which may be randomly distributed (in per- 
sistent sources) or spin-downs in quiescence followed 
by rapid spin-ups during bright transient phases (in 
transient pulsars) or long monotonic spin-up and spin- 
down episodes accompanied by spectral and luminos- 
ity changes (e.g. 4U 1626-67, |Yi fc Vishniac 1999|). 



The AXPs are in many respect also similar to the 
X-ray counterparts of the Soft Gamma-ray Repeaters 
(SGR). The X-ray spectral and timing properties of 
these two type of objects have strong similarities, half 
of the AXPs and SGRs are associated with supernova 
remnants. This has lead to the suggestion that the 
AXPs are also magnetars in which the X-ray emission 



is due to magnetic field decay (Thompson & Dimcan 
"T99^ ). The main difference between these two type of 



objects is the non detection of SGR bursts from the 
AXPs. However, considering the rari ty of the SGR 
activity a mong the established SGRs ( Kouveliotou et 
al. 1996 ), the absence of bursts from AXPs is not 
a serious issue. From a relatively young age of the 
AXP, IE 1841—045 in the sup ernov a remnant Kes 73, 
Gotthelf, Vasisht, & Dotani ( 1999 ) proposed that in 
the evolutionary track, the AXPs are an early quies- 
cent state of the SGRs. Stability of the X-ray emis- 
sion properties (spin-down rate, luminosity, spectral 
shape and pulse shape and fraction) is usually men- 
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tioned as one important aspect of the AXPs though 
one has to compare between observations made with 
different instruments for which the energy band, en- 
ergy resolution and sensitivity are not identical. To 
make a rigorous comparison in the stability of the 
X-ray emission properties we have made new obser- 
vations of two AXPs with the ASCA, four years after 



two previous observations reported by White et al 
(19p)| and Corbet & Mihara ( |1997| ). The aim was to 
critically examine the stability of the X-ray emission 
pattern and more pulse period measurements which 
may provide support to either the accretion powered 
or the magnetar hypothesis for these objects. 

The source 4U 0142-1-61 is close to a long period 



bin hrv pulsar T!X .70146.9+6121 a,nd in the EXOSAT 
observations of this field in which pulsations were 
first discovered, pulsation s from both the sources were 
observed simultaneously ( fsrael, Mereghetti fc Stella 
19c(|). With EXOSAT, the 8.7 s pulsations in this 
source were detected only in the 1.0-4.0 keV range. 
A binary nature of the system was preferred in spite 
of a large X-ray to optical flux ratio and absence of 
pulse arrival time delay. The source was subsequently 
observed with RO S AT a.nd the pulse period was found 



to be very close to the EXOSA T measurement ( Hel- 
lier| 1994|). ASCA observation in 1994 confirmed the 



steady spin-down trend and defined the spectral char- 
acter clearly. A model consisting a 0.4 keV blackbody 
and a power law with a photon index of 3.7 was found 
to describe the spectrum well. From a small radius of 
a few km of the black body emission zone (which prob- 
ably is on the surface of the neutron star) [White ct 



al. (1996) suggested that the black-body component 



is more likely to be due to a spherical accretion rather 
than accretion from a disk. A small pulse fraction and 
energy dependent pulse profile, double peak at low en- 
ergy and single peak at high energy is characteristic 
of this pulsar. From RXTE observations, the upper 
limit on the pulse arrival time delay was determined 
to he 260 ms in the 70 s to 2.5 days range thereby rul- 



ingl out all types of binary companions except white 



dw arf or low mass He main sequence star ( [Wilson ct 
al. [ 1999|). BeppoSAX observations during 1997-98 



confirmed the spectral characteristics, pulse profiles 
and spin-down trend ( Israel et al. 1999h| ). Observa- 
tions spanning 20 years during 1979-1998, with the 
Einstein, EXOSAT, EOS AT, ASCA, BeppoSAX and 
RXTE show an overall constant spin-down trend with 
P = 2.2 X 10-12 s s-i. 

Pulsations in the source IE 1048.1—5937 were dis- 



covered from observations with the Einstein observa- 
tory in 1979 and was confirmed by EXOSAT obser- 
vation in 1985 (Seward, Charles, & Smale 1986). The 
energy spectrum was found to be a power-law type (F 
= 2.26) with low energy absorption (Nh = 1.6 x 10^^ 
atoms cm~2). The relatively harder power-law spec- 
trum and a candidate optical counterpart lead to the 
speculation that this can be a Be star binary. Sev- 
eral CINCA observations established a secular spin- 



few other soft spectrum pulsars ( [Corbet fc Day 1990| ). 
Subsequent i?OS'^T observations however revealed an 
increase in the spin-down rate which is remarkably 



differ ent from other established AXPs ( Mereghett 



1995| ) . The low energy part of the spectrum was first 



accurately me asured with the ASCA in 1994 (Cor- 
bet & Mihara 1997). However, it could not distin- 
guish between power-law (F = 3.34) and a combina- 
tion of power-law (F = 2.0) and black-body (kT = 
0.55 keV). During the ROSAT and the ASCA obser- 
vations, the spin-down rate remained at a higher level 
of 3.3 X 10^11 s s~^. A decrease in intensity by a fac- 
tor of 3 compared to the EXOSAT observation was 
also noticed. BeppoSAX observation in 1997 showed 
that a combined black-body and power-law model fits 
the data well. The size of the black-body emission 
zone was found to be of the order of km^, identical to 
other AXPs ( posterbroek et al. 1998] ). Compared to 
other AXPs, the pulse fraction was found to be much 
larger (- 65%) in IE 1048.1-5937 and has httle en- 
ergy dependence. The power-law component is also 
relatively harder (F = 2 ~ 2.5) compared to the other 
AXPs (F = 3 - 4). RXTE observations in 1996-97 
showed yet another change in the spin down rate, now 
close to that in 1980's. Long RXTE observations es- 
tablished a small upper limit of 60 ms for any pulse 
arrival time delay for orbital period in the range of 
200 s to 1.5 days. Based on this, any binary compan- 
ion other than low mass helium burning star in a face 
on system has been ruled out (Mereghetti, Israel, & 



Stella 1998) 



2. Observations and data analysis 

Both the sources were observed twice with ASCA, 
in 1994 and in 1998 with about 4 years time difference 
between the two observations. ASCA has two Solid- 
state Imaging Spectrometers (SIS) and Gas Imaging 
Spectrometers (CIS) each at the focal plane of four 
identical mirrors of typical photon collecting area 250 
cm^ at 6 keV. The energy resolution is 120 eV and 
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600 eV (FWHM) at 6 keV for the SIS and GIS de- 
tectors respectively. For more details about ASCA 
please refer to [Tanaka, Inouc fc Holt (1994) . Details 
of the 1994 observations are given in Corbet & Mi- 
hara ( |1997| ) and [White et al. (1996)| . In 1998, the 



GIS observations were made in normal PH mode in 
which the time resolution is 64 ms and 500 ms at 
high and medium bit rates respectively. The SIS ob- 
servations were made with one of the CCD chips, and 
has time resolution of 4 s. The standard data se- 
lection criteria of the ASCA guest observer facility, 
that comprises a cut-off rigidity of charged particles 6 
GeV/c, maximum rms deviation from nominal point- 
ing of 0.01 degree, minimum angle from Earth's limb 
10°, satellite outside the South Atlantic Anomaly re- 
gion etc were applied. Data were removed from the 
hot and flickering pixels of the SIS detectors and also 
the charged particle events were removed from the 
GIS detectors based on rise time discrimination. For 
the two observations of 4U 0142-1-61, the source pho- 
tons were extracted from circular regions of radius 6 
arc min and 4 arc min around the source for GIS and 
SIS respectively. In the case of IE 1048.1-5937, the 
source photons were extracted from relatively smaller 
regions of 5 and 3 arc min respectively. For SIS, the 
background spectra were accumulated from the whole 
chip excluding a circular region around the source and 
for GIS it was collected from regions diametrically op- 
posite to the source location in the field of view. 

2.1. Period analysis 

To calculate the pulse periods accurately, barycen- 
tric correction was applied to the arrival time of each 
photon and light curves were extracted from the pair 
of GIS detectors with a time resolution of 0.5 s in the 
energy band of 0.5-10.0 keV. Epoch folding method 
was applied to obtain the pulse periods approximately 
and templates for the pulse profiles were created by 
folding the light curves at the approximate pulse peri- 
ods. Subsequently, the light curves were divided into 
eight segments of equal length and pulse profiles were 
created from each of these segments by applying the 
same epoch and pulsation period. The relative phases 
of the pulses were then evaluated by cross correlating 
the pulse profiles with the respective templates. A lin- 
ear fit to the relative phases with their pulse numbers 
gave the correction necessary to obtain the accurate 
pulse period. The 1998 observations of the two AXPs 
resulted in new measurement of pulse periods at these 
epochs and from the 1994 observations we obtained 




Pulse phase 

Fig. 1. — The background subtracted pulse profiles 
of 4U 0142-1-61 obtained with the GIS in three energy 
bands are plotted for two cycles. A change in the 
pulse profile, from double peaked at low energies to 
single peaked at high energies can be noticed. 



similar pulse periods as reported earlier, with reduced 
uncertainty for 4U 0142-1-61. The pulse periods ob- 
tained for the two sources are given in Table 1. The 
pulse profiles of the two sources in three energy bands 
are shown in Figure 1 and Figure 2. The pulse frac- 
tion, defined as the ratio of the pulsed to total flux, 
was calculated from background subtracted pulse pro- 
files in the 0.5-10.0 keV band. Pulse fractions were 
found to be identical in both the observations, ^ 9% 
and ~ 75% in 4U 0142-^61 and IE 1048.1-5937 re- 
spectively. In 4U 0142-1-61, the pulse profile shows 
energy dependence, from double peaked in low (0.5- 
1.5 keV) and medium (1.5-4.0 keV) energy to single 
peaked in high (4.0-8.0 keV) energy, associated with 
increase in the pulse fraction. In IE 1048.1—5937, on 
the other hand, the pulse profile and pulse fraction are 
almost identical over the ASCA energy range. Light 
curves of the two sources did not show any intensity 
variations at minutes to days time scale. 

2.2. Spectral analysis 

In 4U 0142-1-61, a simple absorbed power-law fit 
shows large residuals at low energy and addition of a 
black-body component results in acceptable fit. For 
the 1998 observation, inclusion of the black-body 
component improved the reduced from 2.88 to 
0.94 and 1.23 to 0.76 for the GIS and SIS respec- 
tively. Similar improvement in fitting was reported 



by White et al. (1996) for the 1994 observation. The 
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Pulse phase 



Fig. 2. — The background subtracted pulse profiles 
of IE 1048.1-5937, similar to Figure 1. The pulse 
fraction of this source is ^75%, highest among the 
AXPs. 




Energy (keV) 



Fig. 3. — The observed SIS and GIS energy spectra 
of 4U 0142+61 shown with histograms for the model 
spectra folded through the responses matrices and the 
residuals. 



photon index of the power-law component obtained 
from the SIS data for the two observations are identi- 
cal, r — 3.3. The photon index from the GIS data is 
slightly different from the SIS value, F = 3.9, but it is 
identical for the two observations. The temperature 
of the black body component obtained from both SIS 
and GIS are identical, 0.39 keV, in both the observa- 
tions and the flux in the two components are also iden- 
tical. The difference in the photon index between the 
SIS and GIS can possibly be attributed to calibration 
uncertainty, because the in-flight spectral calibration 
of the spectrometers are done with sources which have 
relatively harder spectra. The SIS and GIS spectra of 
the 1998 observation are shown in Figure 3 along with 
the best fitted models for the respective detectors and 
the residuals to the model spectra. 

For the other source IE 1048.1-5937, Corbet and 



Mihara (1997) showed that the power-law and black- 
body with power-law models could not be distin- 
guished from spectral analysis of the ASCA data. We 
have found different and deeper minima in the for 
the black-body with power-law model. The improve- 
ment in x^i (^X^ of 41 for 292 degrees of freedom 
and 28 for 261 degrees of freedom for the 1998 and 
1994 observations respectively, both indicating prob- 
ability of chance occurrence less than 10~^) therefore 
favours the inclusion of a black-body component. The 
black-body component has also been detected with 
BeppoSAX ( posterbroek et al. 199^ ). The best fitted 
power-law plus black-body model gives a photon in- 



dex of ^ 3.0, black body temperature of ^ 0.56 keV 
and column density of ~ 1 x 10^^ atoms cm~^. The 
spectral parameters are identical in both the ASCA 
observations and are given in Table 2. Due to relative 
weakness of this source, simultaneous spectral fitting 
was carried out with the GIS and SIS data for the 
1998 observation (Figure 4). However, for the 1994 
observation, the SIS were operated in FAST mode 
and spectra were available only from the two GIS de- 
tectors. 

3. Discussion 

3.1. Period changes and Emission stability 

The pulse period measurements of 4U 0142-1-61 is 
rather scarce except for the last two years (Figure 
5). In spite of the source being very bright, with 
a flux of more than 10"^" erg cm~^ in the 2.0-10.0 
keV band, a pulse fraction of only about 10% has 
restricted the individual pulse period measurements 
accurate to only about AP/P ^ 10~^. A linear fit 
to the pulse period history shows that of the 10 mea- 
surements available, only during the 1994 ASCA ob- 
servation the pulse period measurement was slightly 
different (2.5cr) from the linear trend. If the reported 
errors of all the measurements are taken at their face 
value, the linear fit gives a reduced of 1.7 for 8 de- 
grees of freedom. It therefore can be concluded that 
the recent ASCA measurement together with the pre- 
vious results is consistent with a constant spin-down 
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Energy (keV) 

Fig. 4. — The SIS and GIS energy spectra of 
IE 1048.1-5937 along with the residuals. The his- 
tograms represent the best fitted model folded with 
the response matrices. A simultaneous fitting of the 
SIS and GIS spectra was performed. 



rate. The observations are not yet sufficient to clearly 
identify any significant variation from a constant spin- 
down rate. 

In the source IE 1048.1—5937, departure from a 
linear spin-down is already known ( Mereghetti 199^ ; 
Oosterbroek et al. 1998). In this source, an order of 



magnitude larger spin-down rate and better pulse pe- 
riod measurements (AP/P ~ 10~^), owing to a high 
pulse fraction ('^ 75%), help us to identify three differ- 
ent epochs of spin-down history. Though the method 
adopted for calculating the errors in the pulse pe- 
riod is not known for all the observations and the 
uncertainty level is likely to be nonuniform, a con- 
stant spin-down trend can be ruled out without any 
doubt. A linear fit to the pulse period history with 
the reported errors gives a reduced of 4500 for 11 
degrees of freedom. Including the 1998 ASCA ob- 
seryation (see Figure fi) with the recent Rpppnfi A X 



a.TiH RXTE nbservatinus^ we find that from 1 99fi the 



source has a spin-down rate of (1.67 ± 0.02) x 10 s 
s~^. This is a factor of 2 smaller than the spin-down 
rate of (3.29 ± 0.03) x 10"" s s^^ during the 1994- 
1996 period. The present spin-down rate is closer to 
the value of (1.5 ± 0.5) x 10""'^^ s s""'^, measured dur- 
ing the Einstein, EX OS AT and GINGA observations 
made in the period 1979-1988. The spin-down rate 
is much closer to being constant during these three 
epochs with reduced of 0-8, 7.7, and 36 for 3, 3, 
and 1 degrees of freedom respectively. 



These two sources do not show flux variability on 
time scales from a few minutes to days. In the ASCA 
observations of both the sources separated by 4 years 
we have found that the overall intensity and spectral 
parameters have remarkable stability. A difference 
between the GIS and SIS photon index that has been 
found in 4U 0142+61, is due to calibration uncer- 
tainties. The spectral parameters obtained from the 
1998 GIS and SIS observations are identical to the 
1994 values. The spectral parameters obtained from 
the simultaneous fitting of the GIS and SIS spectra 
are similar to the BeppoSAX values obtained during 
1997-1998. In 4U 0142-^61, the flux history shows a 
rms variation of 15% around the average value (Figure 
5), and multiple measurements with the same instru- 
ment {ASCA and BeppoSAX) gave almost identical 
flux. In IE 1048.1-5937, the over all intensity dur- 
ing the two ASCA observations and one BeppoSAX 
observation in between are within 10% of the average 
value. The 2.0-10.0 keV fluxes during the Einstein 
and ROSAT observations are estimated by extrap- 
olating the measurements in the low energy bands 
of 0.2-4.0 keV and 0.5-2.5 keV respectively, and us- 
ing a rather low photon index of 2.26, obtained by 
EXOSAT. The flux during the GINGA observation 
is estimated by comparing the pulsed fluxes during 
the GINGA and _E'X0S'j4T observations and assuming 
that the pulse fraction remained same. The flux mea- 
surements from the previous observations as shown 
in the bottom panel of Figure 6 are about a factor 3 
higher than the recent measurements with ASCA and 
BeppoSAX. We note that there is some overestima- 
tion in extrapolation of the soft X-ray measurements 
with Einstein and ROSAT dxie to a smaller photon in- 
dex used, and the flux estimates from EXOSAT and 
GINGA could be overestimated due to contribution 
from the nearby bright and variable source 77-Carina, 
which is only 0.4° away and is about 15 times brighter 
than IE 1048.1-5937 dCorcoran et al. 199^; [ishibashi 



et al. 1999) ) . In view of the stability of the flux during 



1994-1998, as obtained from the imaging instruments 
ASCA and BeppoSAX, it is possible that the overall 
intensity of this source does not vary at a few years 
times scale. The spectral parameters from the two 
ASCA observations are also identical and consistent 
with the values obtained with BeppoSAX in between. 
However, we note that in IE 1048.1—5937, the ab- 
sorption column density obtained from the two ASCA 
observations are identical and a factor of 1.5-2 larger 
than the BeppoSAX measurement in between. This 
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Fig. 5. — The pulse period and fLux history of 
4U 0142+61. The straight Une shows the best fit for a 
constant spin-down. The pulse period measurements 
with ASCA are marked with filled circles and the open 
circles are for all the other observations mentioned in 
the text. 



difference is somewhat larger than the known calibra- 
tion difference between the two instruments (Orr et 
al 



1998) 



Some Low Mass X-ray Binaries (LMXB) have also 
been detected at low luminosity levels 10^'^ erg s""'^; 
e.g, Cen X— 4, Aql X— 1, see Tanaka & Shibazaki 1996 



and references therein), understood to be quiescent 
phase of the Soft X-ray Transients (SXT). Usually 
the LMXB sources show both short and long term ir- 
regular intensity variations and many also show quasi- 
periodic oscillations, bursts, dips or orbital modula- 
tion. Even though the AXPs do not show significant 
temporal variations other than the pulsations, which 
is somewhat different from typical characteristics of 
low-luminosity LMXBs, this itself is not a strong ar- 
gument against the AXPs being LMXBs. One LMXB 
which has some properties similar to the AXPs is 4U 
1626—67. The X-ray luminosity, magnitude of spin- 
change rate, pulse period, and flux stability over very 
short to years time scale of this source are similar to 
the AXPs. But, presence of both spin-up and spin- 
down, quasi-periodic oscillations, optically bright ac- 
cretion disk, and a hard X-ray spectrum makes it dif- 
ferent from the AXPs. 



Fig. 6. — The pulse period and flux history of 
IE 1048.1-5937. The lines are used to identify the 
different spin-down epochs. The pulse period mea- 
surements with ASCA are marked with filled circles 
and other observations with open circles. 

3.2. Accretion torque in the common enve- 
lope evolution model 

It has been proposed that the AXPs are recent 
remnants of common-envelope evolution of high-mass 
X-ray binaries (van Paradijs et al. 1995; Ghosh et al. 
1997). In this model, the pulsar is rotating near its 



equilibrium period, close to the Keplerian period at 
the innermost part of the disk. If the pulsar is rotating 
at the equilibrium period, small changes in the mass 
accretion rate cause alternate spin-up and spin-down 
episodes. The over all spin-down is explained with 
the assumption that in the absence of a companion 
star as the source of mass accretion, the mass accre- 
tion rate from the disk decreases slowly on viscous 
time scale. The equilibrium period of the pulsar is in- 
versely related to the mass accretion rate and shows 
secular increase. For IE 1048.1—5937, which has a 
pulse period of 6.5 s and luminosity of 6.3 x 10'^^ erg 
s~^ for a distance of 3 kpc (or in a more favourable 
case, 7 x 10^'* erg s~^ if the source is at a distance 
of 10 kpc), the magne tic field strength inferred for 
an equilibrium rotator ( Frank, King fc Rainc 199^ ) is 
B = 10" (P/3)^/'^ L35^/2 = 6.2 X lO^" (or 2.1 x 10") 
gauss. Here, and in what follows we have assumed 

1.4 M0, 



a neutron star with mass M 



radius R = 

10'' cm, and moment of inertia I = 10^^ gm cm^. A 
pulse period of P implies that if in equilibrium, the 
co-rotation radius or in this case the radius of the in- 
ner disk is rM = (GM)i/3 (p/27r)2/3. The infalling 
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(1. IE 1B41-045 
(3. IE 2259+586 
(3, IRXS 1706-40 
(4. 4U 0142+61 
(5. IE 1046.1-5937; 
(6. SGR 1900+14 
(7. SGR 1806-20 



luminosity or the spectral parameters. 

Among the other AX P s, IRXS J170849. 0-4 00910 
( Sugizaki et al. 1997 ; Israel et al. 1999a ) and 



IE 1841-0 45 (K es 73, [Vasisht k Gotthelf 1997| ; Got- 
thelf et al. 1999| ) have relatively large spin-down rate 



of 2.2 X 10^^ and 4.1 x 10~^^ s respective, while 
the X-ray luminosity is in the range of 10''^ — 10"^^ 
erg s~^. In the common envelope evolution model of 
the AXPs, a faster spin-down compared to the torque 
provided by the accreting matter can be a potential 
problem for these two sources also. For this model 
to be correct for the AXPs, a very fine tuning of the 
mass accretion rate is required. All the sources need 
to have disk accretion rate a tiny fraction larger than 
the propeller regime. Li (1999) have identified sev- 
eral other problems in the context of applicability of 
this model to the AXPs. Most notable is the limited 
lifetime of an accretion disk around a solitary neu- 
tron star compared to the response time scale of the 
neutron star to changes in the accretion torque. One 
possible alternative is that the spin-down is due to 
magnetically driven wind from an accretion disk, pro- 
posed also for SGR 1900-f 14 (|Marsden, Rothschild, & 



0.2 0.5 1 2 5 10 
B (10^"^ gauss) 

Fig. 7. — The X-ray luminosity and pulse fraction 
of the AXPs and SGRs are plotted against the mag- 
netic field strength assuming that these objects are 
magnetars. Two values of the luminosity are plotted 
for IE 1048.1—5937 to show the uncertainty in its dis- 
tance and same has been done for SGR 1900-1-14 to 
sh ow the variability in its X-ray emission. Distance 



of : RXS J170849.0-400910 is taken to be 10 kpc and 



Lingenfelter 1999), but this will require a harder X- 



that for other sources are the best available estimates. 



ray spectrum than what has been observed. Chatter- 



jee, Hernquist fc Narayan (1999) proposed a scenario 



material from the disk carries a positive angular mo- 
mentum of M (GMrivi)^^^. The torque can also be 
expressed in terms of the pulse period P and the X- 

1 /3 

ray luminosity Lx as RLx (P/27rGM) ' . Even if we 
assume that all of the X-ray emission is a result of 
disk accretion, the accretion torque is only 1.1 xlO'^^ 
(or 1.2 xlO'^^) gm cm^ s~^. To achieve the observed 
spin-down rate for a neutron star with moment of in- 
ertia 10^^ gm cm^, the negative torque required to be 
imparted onto the neutron star is iCl — 4.9 x 10'^'^ gm 
cm^ s^^. This is a factor of 450 (or 40) larger than the 
accretion torque, and in the common envelope evolu- 
tion model, a negative dimensionless torque of this 
magnitude is required to spin-down the pulsar at the 
observed rate. In other words, the spin-down rate of 
this source is much larger than what can be achieved 
with disk accretion onto a neutron star with a lumi- 
nosity of less than 10^^ erg s"^ IE 1048.1-5937 does 
not fit in the classical picture of P vs PU^^^ of the 
equilibrium rotators (Ghosh & Lamb 1979). In ad- 
dition, the two ASCA observations made during the 
two epochs which have a factor of 2 different spin- 
down rates, do not show significant difference in the 



in which the AXPs are formed as single neutron stars 
with fossil disks made from fallback material from the 
supcrnovae explosions. In this model, the spin-down 
from an initial period of a few ms to ^6 s is due to 
strong propeller effect at some period when the accre- 
tion rate is very low. But, if accretion is the correct 
phenomenon in the AXPs, as high spin-down goes on 
in presence of substantial accretion, spin-down due 
to wind outflow seems to be more plausible than ac- 
cretion induced angular momentum loss. However, it 
should be remembered that the classical equilibrium 
disk picture assumed here is often found not to be 
the most appropriate des cription for the X-ray pul- 
sars (iBildsten et al. 1997|) . 



3.3. Magnetar model 

In view of a very narrow mass and type allowed for 
any binary companion, and several arguments against 
the common envelope evolution model, the magne- 
tar model seems to be the most likely one for the 
AXPs. If the spin-down is due to magnetic braking, 
dipole field strength of the order of 10^^ gauss is esti- 
mated for these sources. A nearly constant spin-down 
property was thought to favour the magnetar model 
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over a binary scenario. In two AXPs, IE 1841—045 
(Gotthelf et al. |1999| ), and IRXS J1 70849 .0-400910 

(Kaspi, Chakrabarty & Steinberger 1999), there is 
very strong evidence for constant spin-down, whereas 



in IE 1048.1—5937, deviation from a hnear trend 



is clear. Recently, a deviation has also been de- 
tected from SGR 1900+14 (Woods et al. |1999D 



IE 2259+586, the most frequently observed AXP, has 
provided an interesting pulse period history. Obser- 
vations made for about 15 years with many instru- 
ments preceding RXTE showed considerable va riation 
in the spin-down rate ( Baykal fc Swank 1996 ). But, 
the pulse-coherent timing observations with RXTE 
proved it to be othe rwise , at least for a period of last 
2.6 yr (Kaspi et al. |l999|) . 



Two scenarios have been proposed which can ex- 
plain the changing spin-down rate even when the over- 
all braking is d ue to the ultrastrong magnetic field. 
Melatos (1999) showed that for reasonable neutron 
star parameters, a radiative precession effect may take 
pla ce which can give the observed spin-down varia- 



tio Jis with time scale of about 10 vears. It will be 



possible to verify this scenario when more pulse pe- 
riod measurements become available in the next few 
years. Woods et al. ( 199S| ) have identified a possible 
braking glitch in SGR 1900+14 close to the time when 
SGR activity was very strong. But if the spin-down 
variation is related to the SGR activity, similar activ- 
ities should have been observed from IE 1048.1—5937 
and IE 2259+586. We have found that all the gamma 
ray bursts observed with the BATSE for which the es- 
timated positions are within 2a of these AXPs (about 
30 GRBs around each AXP), have strong high energy 
emission unlike the SGR bursts. Heyl & Hcrnquist 
(19 p9)| have proposed that the spin-down variations 
can be explained as glitches (similar to radio pulsars) 
superposed on constant spin-down. But, with recent 
pulse period measurements of IE 1048.1—5937 and 
SGR 1900+14, this will require too many glitches, 
one before almost every observation unless there are 
braking glitches, never observed in radio pulsars. 

In the magnetar model, the X-ray emission is due 
to decay of the magnetic field. The energy gener- 
ated at the core is transported to the crust along 
the magnetic field direction. The black body compo- 
nent of the spectrum is thermal emission from the hot 
spots at the magnetic polar regions and the power- 
law component is part of the thermal emission re- 
processed by the magnetic field and the environment. 
Investigation is required about the expected pulse 



profile and its energy dependence. Time and/or en- 
ergy dependence of the pulse profile as has been ob- 
served in SGR 1806-20 (|Kouveliotou et al. 1998|) and 
SGR 1900+14 ([Hurley et al. 1999| ; [Kouveliotou et al ' 



1999; Murakami et al. 1999) also requires to be ad- 



dressed. A double peaked pulse profile at low energy 
and single peaked profile at high energy are observed 
in 4U 0142+61 (Figure 1). The pulsation is very weak 
in some sources (only 5-10% in SGR 1806—20 and 
4U 0142+61), and 75% in IE 1048.1-5937. If the 
pulsation is due to confinement of the heat in the mag- 
netic polar regions by the magnetic field, a correlation 
between magnetic field strength and pulse fraction 
should be observed. But, it is likely to be smeared by 
the geometric effect of individual sources, i.e. the ori- 
entation of the spin and magnetic axes with respect to 
the line of sight. In the magnetar model, there are two 
mechanisms by which X-rays can be generated. If the 
X-ray emission is powered by decaying magnetic field, 
the luminosity is a very strong function of the mag- 
netic field strength, Lx oc B"* ( [Thompson fc Duncan 
1996). Alternate process of X-ray generation is parti- 



cle acceleration by Alfven waves resulting from small 
scale fracture of the crust, in this case Lx oc B^. The 
X-ray luminosity and pulse fraction of five confirmed 
AXPs and two SGR sources are shown in Figure 7. 
against the magnetic field strength. The later is es- 
timated from pulse period and the overall spin-down 
rate assuming B = 3.2 x 10^^ (PP)^/^ gauss. Though 
there is uncertainty in the luminosity of some sources 
(see the caption of Figure 7), a 2nd or 4th power 
correlation between Lx and B does not seem to be 
present. There is also no correlation between pulse 
fraction and the magnetic field. 

A clustering of the pulse period of 10 sources (7 
AXPs and 3 SGRs) in the 5-12 s range also needs 
to be addressed, when the magnetars are expected to 
be alive in X-ray until they have slowed down to a 
pulse period of about 70 s ( Duncan fc Thompson ). 
The magnetars are expected to be radio quiet due 
t o suppression of pair cre ation at high magnetic field 
( Baring fc Harding 1998[ ), and this seems to be true 
for most sources. 

Changes in luminosity or small changes in the spec- 
tral parameters can be a result of varying activities 
in the core, the heat generated from which is trans- 
ported to the surface along the direction of the mag- 
netic field at time scale of a few years. Variability 
study of AXPs with very high sensitivity may rule 
out magnetar model if significant change in column 
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density is observed. This will indicate the presence 
of accretion disk and wind in the neighbouring area. 
The BeppoSAX and ASCA observations give identical 
column density for 4U 0142+61 but slightly different 
column density in case of IE 1048.1—5937. However, 
multiple observation of the later source with the same 
instrument is found to give identical value indicating 
that the difference between the two instruments can 
also be a systematic effect. 

4. Conclusion 

Using multiple observations of two AXPs with 
the ASCA we have found remarkable stability in 
the intensity, spectral shape and pulse profile over 
a 4 years period. For the source IE 1048.1—5937, 
we have confirmed that similar to other AXPs, the 
spectrum consists of a power-law component and a 
black body component. The spin-down trend of 
4U 0142+61 is consistent with a constant rate whereas 
in IE 1048.1—5937 we have clearly identified three 
different spin-down epochs. We have shown that the 
fast spin-down of some of the AXPs is difficult to 
achieve with disk accretion. Hence the common en- 
velope evolution scenario of AXPs is unlikely to be 
the case, imless the spin-down is due to wind out- 
flow from the disk. In this case also, the stability 
of X-ray emission in spite of varying spin-down rate 
remains unexplained. The significantly different spin- 
down rates of IE 1048.1—5937 in different epochs are 
also difficult to reconcile in the magnetar model, un- 
less precession is at work. But the flux stability and 
lack of variation of the spectral parameters appear to 
favour the magnetar model. 

We thank an anonymous referee for many sugges- 
tions which helped to improve a previous version of 
the manuscript. B. Paul was supported by the Japan 
Society for the Promotion of Science through a fel- 
lowship. 
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Table 1 

The pulse periods from the ASCA observations 



Date of Observation 


Source name 


Pulse period 


(MJD) 




(s) 




4U 0142+61 




49614.1 




8.687873 ± 0.000034 


51046.7 




8.688267 ± 0.000024 




IE 1048.1-5937 




49416.5 




6.446645 ± 0.000001 


51021.1 




6.450815 ± 0.000002 



Table 2 
The spectral parameters 



Obs. date 



GIS 



1994/09/18-19 
SIS 



4U 0142+61 

SIS+GIS 



GIS 



1998/08/21 

sis 



SIS+GIS 



Plioton index 
Power-law norm'' 
BE temp (keV) 
BB norm° 
Reduced / dof 
Observed flux'' 



1.03 ±0.08 
3.9 ±0.1 

0.26 ±0.05 

0.39 ±0.01 
1.4 ±0.2 

0.95 / 691 



0.92 ±0.05 

3.3 ±0.1 
0.12 ±0.03 
0.380 ±0.006 
1.9 ±0.1 
1.11 / 356 
13.0 



1.10 ±0.04 
3.84±0.08 
0.25 ±0.03 
0.382 ±0.007 

1.5 ±0.1 
1.35 / 1030 



1.08 ±0.09 
3.98 ±0.15 
0.30 ±0.06 
0.399 ±0.014 
1.16 ±0.18 
0.76 / 354 



0.97 ±0.08 

3.3 ±0.2 
0.10 ±0.03 
0.384 ±0.005 
1.97 ±0.15 
0.94 / 287 
13.4 



1.17 ±0.04 
3.87 ±0.09 
0.24 ±0.03 
0.378 ±0.009 
1.30 ±0.12 
1.76 /645 



Obs. date 



GIS 



1994/03/02-05 



IE 1048.1-5937 



1998/07/26-27 



SIS+GIS 



Photon index 
Power-law norm'' 
BB temp (keV) 
BB norm'= 
Reduced / dof 
Unabsorbed flux'' 



1.0 ±0.2 
2.9 ±0.3 
5 ±3 10-3 
0.57 ±0.04 
0.05 ±0.02 
0.98 /261 
0.58 



1.21 ±0.24 
3.2 ±0.5 
7 ±4 IQ-^ 
0.56 ±0.06 
0.060 ±0.015 
0.62 / 292 
0.60 



^10^^ atoms cm ^ 

''photons cm" keV" '^ at 1 keV 

•^10-^ photons cm-2 

^10"" ergs cm-2 s"!, 2-10 keV 
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